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ABSTRACT 

We produce and analyse w-band (A ~ 355 nm) luminosity functions for the red and 
blue populations of galaxies using data from the Sloan Digital Sky Survey (SDSS) 
u-band Galaxy Survey (uGS) and Deep Evolutionary Exploratory Probe 2 (DEEP2) 
survey. From a spectroscopic sample of 41 575 SDSS uGS galaxies and 24 561 DEEP2 
galaxies, we produce colour magnitude diagrams and make use of the colour bimodality 
of galaxies to separate red and blue populations. Luminosity functions for eight redshift 
slices in the range 0.01 < z < 1.2 are determined using the 1/V max method and fitted 
with Schechter functions showing that there is significant evolution in M* , with a 
brightening of 1.4 mags for the combined population. The integration of the Schechter 
functions yields the evolution in the w-band luminosity density out to z ~ 1. By 
parametrizing the evolution as p oc (l+z)' 3 , we find that [3 = 1.36±0.2 for the combined 
populations and = 2.09 ± 0.2 for the blue population. By removing the contribution 
of the old stellar population to the w-band luminosity density and correcting for dust 
attenuation, we estimate the evolution in the star formation rate of the Universe to 
be (3sfr = 2.5 ± 0.3. Discrepancies between our result and higher evolution rates 
measured using the infrared and far-UV can be reconciled by considering possibilities 
such as an underestimated dust correction at high redshifts or evolution in the stellar 
initial mass function. 

Key words: surveys - galaxies: evolution - galaxies: fundamental parameters - 
galaxies: luminosity function, mass function - ultraviolet: galaxies. 



1 INTRODUCTION 

The evolution in the star formation rate (SFR) of the 
Universe over the course of cosmic history is a measure 
of great importance that has many applications in var- 
ious fields of astrophysics. In addition to the obvious 
applications to the build-up of stellar mass in galaxies, 
other examples i nclude constraining the stel l ar initial mass 
funct ion (IMF) l|Baldrv fc Glazebrookl I200&1 ; IWilkins et ail 
120081). testing models of ch emical evolution in the Universe 
(|Calura fc Matteuccil l2004|h constraints on the nature of 



both type l a supernovae (Strolgcr et al 
ray bursts l|Daigne et all 120061 : iLil 120081 ) 



20041 ) and gamma- 
and comparisons 
with global AGN activi t y and the growth of black holes 
jFranceschini et alj[l999l ; Isomerville et alll2008h . 

SFRs of galaxies have been measured by a vari- 
ety of methods, at different wavelengths. These include: 
SFRs determined from the far- and mid-infrared produced 
by the re-emission of ultraviolet (UV) p hotons by dust 
grains, revealing hidden star formation (|Le Floc'h et al.l 
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120051 : iPerez-Gonzalez et al.ll2005h ; SFRs measured from ra- 
dio emission pr oduced by relativistic electrons in super- 
nova remnants (|Condonl I199I ; ISevmour et ai1l200Sft ; Neb- 
ular emission lines such as H a , and [Oil], produced 
by the rec ombination of ionized gas surrounding hot young 
OB stars (iKennicutt fc Kentj|l983l; iMqustakas et al.l l2006l : 
I Cooper et all 120081 : iJames et all l2008al h~ and X-ray emis- 
sion produced by X-ray binary systems in late-type galaxies 
|Norman et al.ll2004l : iLehmer et al]|2008h . 

More directly, SFRs have been measured from the UV 
luminosities of galaxies produced from sh ort lived OB stars 
(lifet imes ~ 10 7 yr ) . Earl y studies such as lLillv et al l (|l996l ) 
and iMadau et all |l996) used this to establish that the 
volume-averaged SFR of the Universe has been in decline 
since z ~ 1 to the present day. Although later studies 
|Madau et al.lll99Sl : ISteidel et al]|l999l : I Wilson et all l2002h 
measured the decline with improved precision, it is still 
poorly constrained and its cause is still a major problem 
in the field of observational cosmology. 

As an alternative to the far UV luminosity, the ti-band 
(A ~ 355 nm) luminosity of a galaxy, predominantly emitted 
by young stars of ages < 1 Gyr, has been shown to be a 
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reasonable tracer of star formation (jHopkins et al.l l2003h . 
Using ti-band has two main advantages over the far UV, 
because i) more data are available to i mprove the st a tistics 
and ii) it is less affected by extinction. iDriver et al.l (|2008T ) 
find that at z — 0, the fraction of photons to escape through 
the dust from the local galaxy population is ~ 45 per cent 
in the it-band as opposed to ~ 25 per cent at 200-nm. One 
disadvantage of the it-band as a tracer of star formation 
is the contribution from old stellar populations, which is 
dominant in early-type galaxies and small in star forming 
galaxies. This has to be accounted for. 

So far it-ban d field luminosity fu nctions (LFs) have 
been produced bv lBlanton et all l|2003bl ) for Data Release 2 
(D R2) of the Sloan Dig ital Sky Survey (SDSS) main sample, 
bv lBaldrv et all l|2005l ) for t he SDSS n-band Galaxy Survey 
(wGS) and more recently bv lMontero-Dorta fc Pradal (120081 ) 
using DR6, probing volumes out toz~ 0.2. In this paper we 
extend the analysis to higher redshift, by combining samples 
of galaxies taken from the SDSS uGS and the Deep Evolu- 
tionary Exploratory Probe 2 (DEEP2) survey, to produce 
LFs in the range 0.01 < z < 1.2, which represents ~ 8 Gyr 
of cosmic history from when the Universe was ~ 5 Gyr old. 
Selected by TJ-band, the DEEP2 survey is ideally suited for 
comparison to the SDSS uGS as 7?-band magnitudes at z ~ 1 
correspond to rest it-band magnitudes. We go on to produce 
colour magnitude diagrams (CMDs) fro m the data and make 
use of the observe d colour bimodality (|Strateva et al.ll200ll ; 
Baldry et al. 20041) to separate the narrow band of red, pas- 
sive galaxies and blue star forming populations of galaxies. 
Integration of the LFs produced from these separate popula- 
tions allows the luminosity density, or total output of tt-band 
light in a given volume to be determined. 

The structure of this paper is as follows: Section 2 de- 
scribes the DEEP2 Survey and SDSS uGS, and the data re- 
duction employed to obtain our final galaxy samples. In sec- 
tion 3 the methods required to calculate the luminosity func- 
tions including completeness corrections and K corrections 
are described, and we present rest-frame CMDs displaying 
how the red and blue populations are separated. In section 
4 we present the luminosity functions, fitted Schechter pa- 
rameters and luminosity densities. Our results are then com- 
pared with the observations at different wavelengths and we 
discuss our estimate of the evolution of SFR of the Universe 
in section 5. Section 6 contains a summary of our conclu- 
sions. Throughout this paper we adopt a cosmology with 
parameters (fl m ,flA)o = (0.3,0.7) and Ho = 70 km s _1 
Mpc -1 . Magnitu des are corrected for G alactic extinction us- 
ing dust maps of ISchlegel et~ail l|l99Sl) . 



2 DATA 
2.1 DEEP2 

The Deep Evolutionary Exploratory Probe 2 (DEEP2) sur- 
vey l|Davis et alj l2003h is a redshift survey primarily de- 
signed to probe the properties of galaxies at redshift z ~ 
1. When fully complete the spectra of around 50 000 galax- 
ies with redshifts z < 1.4 will have been measured over an 
area of sky totalling 3.5 deg 2 . 

In this paper we use data from the third data re- 
lease, DR3, which comprises two catalogues: photometric 



and spectroscopic. The photometric catalogue contains the 
B,R and I magnitudes of 716 465 objects (including stars 
and duplicat e objects) obtained w ith the 12A; x 8fc mo- 
saic camera (|Cuillandre et al.ll200ll ) on the 3.6 m Canada- 
France-Hawaii- Telescope. The spectroscopic catalogue con- 
tains redshifts and spect ra of 46 337 object s measured by the 
DEIMOS spectrograph (|Faber et alj|2003j ) with wavelength 
range 600 — 920 nm and resolution A/ A A ~ 5000 on the 
10m Keck 2 telescope. Galaxies selected for redshift mea- 
surements are limited to apparent AB magnitudes in the 
range 18.5 < R < 24.1. 

The survey area is separated into 4 regions of the sky 
and is made up of a 0.5 sq. deg. field known as the Extended 
Groth Strip (EGS) and three 1.0 sq. deg. fields designated 
fields 2, 3 and 4. Galaxies in fields 2, 3 and 4 are preselected 
by their B, R and I magnitudes using colour cuts, to have 
estimated redshifts z > 0.7. No such preselection was carried 
out for the EGS which means that there is an approximately 
equal number of galaxies below and above z = 0.7 in this 
region. 

2.2 DEEP2 sample reduction 

Before determining the luminosity functions, stars and du- 
plicate objects had to be removed from the DEEP2 cata- 
logues and it was also useful to cut the area of the pho- 
tometric catalogue, due to the spectroscopic catalogue not 
yet covering the same area as the photometric catalogue. 
The methods used to reduce the catalogues to obtain a final 
sample of galaxies used in this paper are described in this 
section. 

The first reduction of the photometric catalogue was 
to cut down the areas of the photometric survey in each 
field, so that the areas of the photometric and incomplete 
spectroscopic catalogues are the same. 177 680 objects in the 
photometric catalogue were cut, reducing it to 538 785. 

The second procedure was to discard galaxies with mag- 
nitudes outside the range 18.5 < R < 24.1, to give the 
spectroscopic and photometric catalogues equal magnitude 
limits. This removed 313 362 objects from the photometric 
catalogue and 58 objects from the spectroscopic catalogue. 
After matching up the spectroscopic object numbers with 
those in the photometric catalogue, 6 objects were found 
to have no corresponding object number in the photometric 
catalogue. These 6 objects were removed from the spectro- 
scopic catalogue. After this reduction the photometric and 
spectroscopic catalogues contained 225 423 and 46 279 ob- 
jects respectively. 

The third reduction was to remove stars from the cata- 
logue. Objects in the photometric catalogue are designated a 
quantity, PGAL (|Coil et al.ll2004h , that gives the probability 
of an observed object being a galaxy or not. This probability 
is determined from star-galaxy separations calculated from 
the magnitudes, sizes and colours of the objects. Resolved 
sources i.e. those which are definitely galaxies are given a 
PGAL value of 3 whereas unresolved sources are given a 
probability of being a galaxy between and 1. In this inves- 
ti gation it was d e cided to use objects with PGAL > 0.2 as 
in lwillmer et all < |2006h . 

The final reduction was to remove galaxies in the cat- 
alogues that had been observed multiple times and desig- 
nated different object numbers. These duplicated observa- 
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Table 1. Properties of the DEEP2 reduced survey sample 



Field 


Reduced Area (sq. deg.) 


N p 


N z 


EGS 


0.45 


36319 


13440 


Field 2 


0.68 


34976 


9559 


Field 3 


1.00 


52390 


11904 


Field 4 


0.94 


42914 


10747 


Total 


3.07 


166559 


45650 



tions for objects were identified by running a procedure 
to find and match objects within 0.5 arc-seconds of each 
other. The matched objects with the best redshift quality 
and smallest error in R mag, were used to produce the lu- 
minosity functions. The duplicate objects with poorer data 
were then removed from the photometric and spectroscopic 
catalogues. These two further processes reduced the photo- 
metric catalogue by 58 824 to 166 599 galaxies. 629 objects in 
the spectroscopic catalogue which had object numbers cor- 
responding to galaxies removed from the photometric cata- 
logue, were also removed, leaving 45 650 objects in the spec- 
troscopic catalogue. 

The properties of the reduced fields can be seen in Table 
[T] Np and N z are the numbers of galaxies with photome- 
try and spectroscopy. Our primary sample consists of 24 561 
galaxies with redshifts in the range 0.4 < z < 1.2. 



2.3 The SDSS u-band Galaxy Survey 

The Sloan Digital Sky Survey is only brie fly describe d here , 
for greater detail the reader is referred to lYork et al.l (|2000l ) 
and subsequent data re lease papers such as the EDR paper 
|Stoughton et al.l 1200 The SDSS is both a photometric 
and spectroscopic survey which makes use of a dedicated 
2.5-m telescope situated on Apache Point, New Mexico. Pho- 
tometry in 5 broadb and filters, ugriz, is conducted using a 
mosaic CCD came ra dGunn et al.l 19981 ) and calibrated with 
a 0.5-m telescope l|Hogg et alj 2001 ). Spectra of selected ob- 
jects are measured using a 640-fibre fed spectrograph with 
wavelength range 380-920 nm and resolution A/ A A ~ 1800. 

The sample of galaxies used in this paper is 
taken from 'Stripe 82' of the SDSS 'Southern Survey' 
ijAdelman-McCarthv et~aD l2006h . This consists of a 2.52° 
wide strip across the Southern Galactic Pole (SGP) cen- 
tred along the celestial equator from RA —50.8° to 58.6°. 
Multiple imaging passes of the stripe allow co-added it-band 
Petrosian fluxes to be determined. This gives an increase in 
the S/N ratio by a factor of about 3 over single epoch obser- 
vations. Ex tra spectroscopic o bservations of galaxies (fully 
described in lBaldrv et alj|2005l ) include a low-z r pelro < 19.5 
sample and a u ae i ec t < 20 sampleQ The multiple scans, and 
extra spectra, allow a catalogue of objects to be produced 
with iipetro < 20.5 based on the co-added data. The result is 
a photometric catalogue of 75 087 objects, with magnitudes 



1 "seiect = u model - r mode i + r petro and can be regarded as 
pseudo u-band Petrosian magnitude. This selection magnitude 
was used because the spectroscopic observations were targeted 
on single epoch data, where u pe tro is less reliable at these faint 
limits. 



in the range 14.5 < u < 20.5 and spectroscopy for 48 640 
objects in an area of 275 deg 2 . Our primary sample consists 
of 41 615 galaxies with redshifts in the range 0.01 < z < 0.2. 



3 METHODS 

3.1 Completeness corrections 

Completeness corrections are required to correct the lumi- 
nosity functions for the number of galaxies that are not tar- 
geted for spectroscopic observations. For both surveys, em- 
pirical completeness values were determined by binning in 
observed quantities. For the SDSS uGS sample of galaxies 
it is assumed that all the redshifts are correct but for the 
DEEP2 data we additionally take into account the redshift 
success rate, which is explained below. 



3.2 DEEP2 completeness 

Redshifts in the DEEP2 survey were first determined by 
an automated pipeline (Cooper et al. 2009, in preparation) 
which does all the image reduction and extraction of spectra 
which are then validated by eye. In this process the redshifts 
were given a quality assessment value (Q) from 1 to 4. Q=l 
indicates the redshift determined to be poor due to a lack 
of features in the spectra, whereas Q=4 means the redshift 
has two or more identified features (the [Oil] A = 372.7 nm 
doublet is considered to be two features) and the determined 
redshift is rock solid. 

In this paper we only use galaxies with Q ^ 3. To cal- 
culate the redshift success, galaxies are divided into bins of 
R band apparent magnitude from 18.5 < R < 24.1. Brighter 
galaxies with 18.5 < R < 20.0 are divided into 0.5 mag bins 
whereas galaxies with 20.0 < R < 24.1 are divided into 0.2 
mag bins. The redshift success is then simply the number of 
galaxies with Q ^ 3 divided by the total number of galaxies 
per magnitude bin. The redshift success as a function of R 
band apparent magnitude is shown in Fig. [T] 

On inspection of Fig. [1] it is clear that the redshift suc- 
cess drops unexpectedly for the brightest galaxies. After cor- 
respondence with Mike Cooper of the DEEP2 collaboration, 
we have found this is mainly due to some low redshift galax- 
ies that have been incorrectly designated as having too poor 
quality spectra (to measure a redshift). The DEEP2 team 
acknowledge this is something that will be corrected for the 
next data release (DR4). Our completeness corrections as- 
sume that the missed redshifts, both bright and faint, have 
a similar redshift distribution as the obtained redshifts. The 
overall conclusions of this paper are unlikely to be signifi- 
cantly affected by the validity of this assumption because 
the completeness is > 0.7 up to R ~ 23. 

A targeting completeness correction accounts for the 
galaxies that are not targeted for spectroscopy. This is 
calculated by dividing galaxies in both photometric and 
spectroscopic catalogues into colour-colour bins. Firstly the 
galaxies are divided into 39 R — I colour bins between 
—3.0 < R — I < 4.0. Galaxies are divided in the follow- 
ing bins [< -0.03,-0.03 -» 1.01 (binsize 0.04), 1.01 -> 1.56 
(0.05), > 1.56]. Then the galaxies are subdivided into 61 
B — R colour bins. All galaxies lie within the range —1.6 < 
B - R < 7.0 and are divided into bins [< -1.0, -1.0 -> 4.9 
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Figure 1. Redshift success fraction as a function of apparent R- 
band magnitude for DEEP2 DR3. The redshift success fraction 
is simply the number of galaxies with redshift quality Q 3 
divided by the total number of galaxies, within an apparent R- 
band magnitude bin. Error bars indicate Poisson errors. 



(binsize 0.1), > 4.9]. This is done with the intention of ob- 
taining a few tens of galaxies in each colour-colour bin. 

Targeting completeness is simply the number of galaxies 
from the spectroscopic catalogue divided by the number of 
galaxies in the photometric catalogue in each colour-colour 
bin. Fig. [2] shows the distribution of colours for all 166 599 
DEEP2 galaxies in the photometric catalogue represented as 
contours along with the targeting completeness as a function 
of R — I and B — R. From this plot the colour preselection 
used to obtain galaxies at z > 0.7 can be clearly seen as 
the region of moderate completeness (green). In this paper 
we choose not to treat galaxies within the Extended Groth 
Strip (EGS) separately from the other DEEP2 fields and an 
effect of this is the region of galaxies with low completeness 
(red) also seen in the figure. 

The overall completeness for a single galaxy is simply 
the product of the redshift success rate and targeting com- 
pleteness. 



3.3 SDSS uGS completeness 

Completenes s corrections f or the SDSS uGS sample are fully 
described in iBaldrv et all (|2005l ). These are determined by 
dividing the sample into bins using the star-galaxy separa- 
tion, A sg (= rpsF — r mo dei, the difference between r-band 
point spread function and 'model' magnitudes), along with 
u,u — g and j-r as variables. Firstly the objects are divided 
into weakly and strongly resolved sources at the dividing 
line of A sg = 0.3. For u > 19.7 strongly resolved sources 
are divided into 0.1 mag bins and weakly resolved sources 
are divided into 0.2 mag bins. These magnitude bins are 
then further divided on the basis of u — g and g — r colour 
bins so that there are around 50 galaxies in each bin. The 
completeness is simply the fraction of photometric objects 
that have spectroscopy. The completeness as a function of 
g — r and u — g colours for galaxies with 19.7 < u < 20.0 is 
presented in Fig. [3] For a more detailed figure showing the 



completeness as a function of A s 



and 



Figure 2. Targeting completeness of DEEP2 galaxies as a func- 
tion of R — I and B — R colours. The targeting completeness 
corrects the LFs for galaxies not selected for spectroscopic obser- 
vations. It is the number of galaxies in the spectroscopic catalogue 
divided by the number of galaxies in the photometric catalogue, 
within a colour-colour bin. All four DEEP2 fields are treated the 
same way regarding targeting completeness. A consequence of this 
is two distinct regions; a region of moderate completeness consist- 
ing of galaxies from fields 2, 3 and 4, having preselected colours to 
have redshifts greater than 0.7 and a region of low completeness 
made up of galaxies from the Extended Groth Strip. The distribu- 
tion of colours for 166 599 DEEP2 galaxies from the photometric 
catalogue are represented by logarithmically-spaced contours. 



1.2 
1.0 



0.2 
0.0 



jmpletenesK 
0.0 - 0.1 
0.1 - 0.2 



Figure 3. Sample completeness as a function of g — r and u — g 
colours for the SDSS «GS galaxies with 19.7 < u < 2 0.0. For com- 



lctcncss as a function of A s 



plctcm 
|2005h 



- g and g — r see IBaldrv et al.l 



Fig. 6 of IBaldrv et al.l ((2005) . Although the completeness 



can be low (~ 0.1-0.2) for the faintest blue galaxies with 
20.3 < u < 20.5, all parts of galaxy colour-colour space over 
the full magnitude range are sampled. 



3.4 K corrections and colour magnitude diagrams 

To convert the observed apparent magnitudes of the galax- 
ies to absolute rest-frame magnitudes we ma ke use of the K 
correct ion program Kcorrect (version 4.1.4) of lBlanton et al.l 
(|2003al l. For the DEEP2 sample, B, R and I magnitudes (ef- 
fective wavelengths « 439, 660, 813.5 nm) are transformed 
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Figure 4. Colour magnitude diagrams for the four SDSS redshift slices, in the range 0.01 < z < 0.2. The number densities are weighted 
by 1/ (VmaxCi) and are represented by logarithmic contours. Little evolution is seen in the colour corresponding to the region between 
red and blue galaxies for all redshift slices. The solid line represents the straight line given by eq.l, that is used to divide the red and 
blue galaxy populations. The dashed lines represent the magnitude limit of the redshift slice. 




-21 -20 -19 -18 -17 -16 -21 -20 -19 —IB -17 -16 



Figure 5. CMD contour plots for four DEEP2 redshift slices in the range 0.4 < z < 1.2. Number densities are represented as logarithmic 
contours. The solid line represents the straight line given by eq.l that is used to divide the red and blue galaxy populations. The dashed 
lines represent the magnitude limit of the redshift slice. 



© 2009 RAS, MNRAS OOO.HTfTil 



6 M. Prescott, I. K. Baldry and P. A. James 
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Figure 6. Evolution of the «-band luminosity function for the 
combined red and blue populations of galaxies. LFs are produced 
for four redshift slices between 0.01 < z < 0.2 for the SDSS uGS 
data and four redshift slices between 0.4 < z < 1.2 for DEEP2, 
using the non-parametric l/Vmax method. Galaxies are binned 
in 0.2 mag bins. 



into rest frame SDSS u and g magnitudes (effective wave- 
lengths 355 and 467 nm). 

For the SDSS uGS sample, the observed u and g mag- 
nitudes are converted into the rest frame. This process en- 
abled the production of rest frame u-g colour magnitude di- 
agrams (CMDs) for the galaxies in 8 different redshift slices. 
Fig. 0] shows CMDs produced from SDSS data and Fig. [5] 
shows CMDs from DEEP2. These are represented as contour 
plots which are weighted by l/(V ma xCi). Here Vmax is the 
maximum comoving volume, within which a galaxy could 
lie within a redshift slice a nd within the magnitude limits of 
the survey (|Schmidtll 968) and d is the completeness. These 
plots show the build up in red galaxies over time and that the 
col our bimod a lity o f galaxies extend s out t o z ~ 1, as found 
bvlBell et al.l il2004ft. IWillmer et all il2006T). Icirasuolo et all 
|2007D . lFranzetti et alj (|2007l ) and lCowie fc Barged (|2008T ). 

On inspection of the DEEP2 and SDSS tfGS contour 
plots, we find although both blue and red populations are 
brightening in it-band magnitude as redshift increases, there 
is no obvious evolution in u — g colours for the natural di- 
viding line. This l ack o f evolution has also been observed by 
ICowie fc Barged (120081) in AB340 0-AB4500 colours out to 
z = 1.5 and lWillmer et all (120061 ) in U - B colours. Given 
this lack of change, red and blue galaxies at all redshifts can 
be separated by a straight line with the equation: 



-0.033M„ + 0.75 



(1) 



Changing this dividing line will have little effect on the over- 
all results, as long as it remains in the gap between blue and 
red populations. 



4 RESULTS 

4.1 Luminosity functions 

The luminosity function (LF) is defined as the number of 
galaxies per unit volume per unit luminosity interval dl, or 
alternatively magnitude interval, AM. In this paper we cal- 



Figure 7. Evolution of the «-band luminosity function for the 
red population of galaxies. 0.2 mag bins are used to produce all 
the LFs apart from the 0.4 < z < 0.6 and 0.6 < z < 0.8 DEEP2 
redshift slices where 0.4 mag binning is used. 
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SDSS 0.10 < z < 0.15 
SDSS 0.15 < z < 0.20 
DEEP2 0.4 < 7 < 0.6 
DEEP2 0.6 < z < 0.8 
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DEEP2 1.0 < z < 1.2 
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Figure 8. Evolution of the u-band luminosity function for the 
blue population of galaxies. 0.2 mag binning is used throughout. 



culate LFs by the 1/Vmux method, where Vmax is defined 
above. This method is non-parametric and does not require 
any assumptions on the shape of the LF. It can be biased 
however if there are strong density inhomogeneities in the 
field. The luminosity function is estimated as follows with 



the number of galaxies per Mpc 
bin given by: 



1 ^ 



AM ^ C,K 



per absolute magnitude 



(2) 



where d is the completeness for each galaxy described in 
the previous section. 

LFs were produced for the red, blue and combined pop- 
ulations of galaxies, for the same 8 redshift slices as the 
CMDs above. For the SDSS uGS data we produced LFs for 
4 redshift slices over the range 0.01 < z < 0.2. For the 
DEEP2 data we produced LFs for 4 redshift slices between 
0.4 < z < 1.2. The data were binned in u-band mag bins 
of width 0.2 mags, apart from the DEEP2 red population, 
0.4 < z < 0.6 and 0.6 < z < 0.8 slices, where 0.4 mag 
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binning was used. The properties of the redshift slices in- 
cluding the numbers of galaxies in the different populations 
and volumes can be seen in Table [2] Figs. [6] [7] and [8] show 
the binned LFs for the 8 redshift slices plotted together for 
the combined, red and blue populations of galaxies. The LFs 
of all the different populations show a shift to higher lumi- 
nosities at greater redshift. 

By binning the galaxies into 8 redshift slices it ensures 
that there is reasonable resolution of the evolution (~ 10 
per cent in luminosity per SDSS slice and ~ 20 per cent in 
the DEEP2 slices) . It also ensures there are sufficient galax- 
ies to produce LFs with good signal to noise after dividing 
into the red and blue populations. This is something par- 
ticularly important for the 0.4-0.6 DEEP2 slice which has 
fewer galaxies. 

For a n analysis of the evolution of the LFs, Schechter 
functions ( Sch echter! Il976h were fitted to the binned data. 
The Schechter function is parametrized as: 

(3) 

where: L* (M*) is the characteristic luminosity (magnitude) 
of a galaxy at the 'knee' of the function; (f>* is the number 
density of galaxies around L* and a is the faint-end slope of 
the function. 

Schechter functions are fitted to the LFs using a least- 
squares fitting method to find the parameters <j>, M* and 
q. In order to quantify the evolution we choose to fix the 
faint-end slope of the LFs for the different populations on 
the assumption that it remains constant with redshift. This 
is due to the faint-end slope being poorly sampled in higher 
redshift slices. The value we take as the faint-end slope is the 
one that minimises the total \ 2 over & U redshift slices. After 
doing this we find best fitting faint-end slopes of a = —1.3 
for the blue population, a — —0.35 for the red population 
and a — —1.0 for the combined LFs. 

Using a slope of a = —1.0 for th e combined population 
LFs is consistent with the results of iBlanton et al.l l|2003bl ) 
who used a sample of 22 020 galaxies with u < 18.4 and 
in the redshift range 0.02 < z < 0.14 from SDSS DR2, to 
obtain a slope of a = -0.92 ± 0.07. iBaldrv et ail (|2005h 
determined a = —1.05 ± .08 fo r z < 0.06 and more recently 
iMontero-Dorta fe Pradal |2008l ) used a much larger sample 
of 192 068 galaxies from SDSS DR6, with 16.45 < u < 19.0 
and a redshift range 0.02 < z < 0.19, to obtain a slope of 
a = -1.01 ±0.03. 

Using photometr ic redshifts from the COMBO- 17 Sur- 
vey |\VbIfetai] (|2003h . produce near-UV (280- nm) LFs for 
25 000 galaxies in the range 0.2 < z < 1.2, as a function 
of SED type. This study finds no evolution in the faint-end 
slope for 4 different SED types. For reddest SED type 1 
galaxies, the faint-end slope is found to be a — —0.5 ± 0.2 
and for the the bluest, SED type 4 galaxies a slope of 
a = — 1.5 ± 0.06 i s foun d. 

Illbert et all (|2005T ) using first epoch data from th e 
VIMOS-VLT Deep Survey (VVDS) l|Le Fevre et al.ll2004h . 
measure an evolution in the faint end-slope of the [/-band LF 
and find that it steepens from a — — 1.05±0.05 in the range 
0.05 < z < 0.2 to a = -1.44±0.2 in the range 0.8 < z < 1.0. 
In our study we do not observe this steepening, but in fact 
the opposite, with z = 0.4-0.8 DEEP2 redshift slices having 
shallower slopes than the SDSS slices (if fitted separately). 



While allowing a to vary provides a better fit to the 
data, it should be noted that the galaxies with the faintest 
magnitudes within each slice have the largest volume and 
completeness corrections. In particular, the faint DEEP2 
data have a low redshift success rate and the assumption 
that the missed redshifts have the same distribution is a 
significant source of uncertainty (cf. the 'average' model of 
iFaber et al]|2007f ). The focus of this paper is on extracting 
reliable luminosity densities, and evolution in M* is more 
robustly determined using fixed a given the well-known de- 
generacy in the M*-a plane. 

The Schechter function fits for the red, blue and com- 
bined populations can be seen in Figs. [9] and 1101 Errors on 
the number density used in the fitting process and shown in 
the figures are Poisson errors added with an error of 5 per 
cent of the number density, in quadrature, to take into ac- 
count systematic uncertainties in the binning process. Dot- 
ted lines in the figures represent the regions outside the mag- 
nitude ranges used to fit the Schechter function. It is notable 
that while a constant value of a is suitable to describe the 
evolution of the blue and combined population LFs, there is 
a deficit of faint red galaxies below the Schechter fit in the 
higher-redshift DEEP2 data (cf. evoluti on of the field dwarf- 
to-giant ratio, iGilbank fc Baloghl [2008 ) . The main focus of 
this paper, however, is on the blue sequence and its rela- 
tionship to cosmic SFR evolution. Tables [3] [4] and [5] display 
how the fitted Schechter parameters evolve along with the 
luminosity density (described in the next subsection) from 
an integration of the LF (pfit) and from a summation of 
the luminosities of all the galaxies (psum)- 



4.2 Luminosity density 

After fitting Schechter functions to the binned LFs, these 
can be used to estimate the total comoving luminosity den- 
sity (LD), with no extinction correction, in each of the red- 
shift slices for the red and blue populations. We assume that 
the fitted Schechter function is valid outside the fitted mag- 
nitude range. The LD is then given in magnitudes per Mpc 3 
by 

j = M* - 2.51og[(07Mpc- 3 )r / (a + 2)] (4) 

where Tf is the gamma function, and in linear units by 

p L = io' 34 - 1 -^)/ 2 - 5 W Hz^Mpc- 3 (5) 

Here j is in AB mag Mpc -3 . 

Fig. [11] shows the evolution in the integrated Schechter 
fit LD for the different populations. For the the blue and 
combined populations, the luminosity density has decreased 
steadily from redshift 1.2 to the present day. The red pop- 
ulation luminosity density appears t o be almost cons t ant, a 
result also found bv lBell et al"Ii|2004l ) and lFaber et al.1 (|2007h 
who conclude this can be explained by passive fading being 
compensated by a build up in the stellar mass density. Note 
also that the fit to the evolution for the blue population LD 
has a lower y 2 than the fit for the red population. This is 
not surprising because the cosm ic variance is larger for the 
more clustered red population (|Somerville et alTl2004l ). In 
fact if the larger error bars from psum (Table 3| are used 
when fitting, the best fit shows a marginally rising LD with 
z rather than declining as in the figure. 
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Table 2. Spectroscopic sample properties 



Redshift Slice 


Sample 


Nall 


N RED 


Nblue 


Volume (10 6 Mpc 3 ) 


0.01-0.05 


SDSS 


6174 


1708 


4466 


0.263 


0.05-0.10 


SDSS 


14002 


4414 


9588 


1.781 


0.10-0.15 


SDSS 


13559 


4302 


9297 


4.614 


0.15-0.20 


SDSS 


7840 


1899 


5941 


8.552 


0.40-0.60 


DEEP2 


1955 


316 


1639 


2.184 a 


0.60-0.80 


DEEP2 


6914 


1296 


5618 


3.404 


0.80-1.00 


DEEP2 


9818 


1363 


8455 


4.488 


1.00-1.20 


DEEP2 


5874 


616 


5258 


5.374 



"This volume applies to all 4 DEEP2 fields, although most galaxies in this redshift slice are contained within the EGS. 
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Figure 9. Luminosity functions of the four SDSS itGS redshift slices, for the red, blue and combined populations of galaxies, produced 
using the 1/V m ax method. The dotted lines show an extrapolation of the Schechter function beyond the fitted magnitude range. The 
LFs for the red and blue populations are scaled down by a factor of 10. Errors on the data points are Poisson errors added in quadrature 
to an error of 5 per cent of the number density, to take into account systematic uncertainties in the binning process, a is fixed at —0.35, 
— 1.3 and —1.0 for the red, blue and combined populations respectively. 



As well as calculating the LD from integration of the 
Schechter fit we also estimate the LD from the summation 
of the galaxies' it-band luminosity divided by l/(V m axCi), 
in the different redshift slices. 



RSUM = ( y ) (6) 

i v / 

The values of this summation for the combined, blue and 
red populations can be found in Tables [3j [4] and [5] 
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Figure 10. Luminosity functions for the four DEEP2 redshift slices, for the red, blue and combined populations of galaxies. For details 
see Fig. [9] Note the deficit of faint red galaxies below the fit for three redshift slices because a is constrained by the lowest two SDSS 
slices. 



Table 3. Schcchtcr function parameters for all galaxies with an assumed a ■■ 



-1.0. 



< z > 


N GAL 


M* 


<t>* (10~ 4 gal Mpc" 3 ) 


logpFiT (W Hz" 1 Mpc -3 ) 


log PSUM (W Hz -1 Mpc -3 ) 


0.030 


6174 


-18.89 ± 0.03 


103.68 ± 2.55 


19.21 ± 0.01 


19.23 ± 0.01 


0.075 


14002 


-19.12 ± 0.02 


81.60 ± 1.79 


19.19 ± 0.01 


19.17 ± 0.05 


0.125 


13599 


-19.27 ± 0.02 


80.57 ± 2.55 


19.25 ± 0.01 


19.10 ± 0.03 


0.175 


7840 


-19.44 ± 0.02 


73.96 ± 3.78 


19.28 ± 0.02 


19.00 ± 0.05 


0.500 


1955 


-19.80 ± 0.05 


69.71 ± 2.74 


19.41 ± 0.02 


19.42 ± 0.05 


0.700 


6914 


-20.07 ± 0.03 


63.71 ± 1.99 


19.46 ± 0.01 


19.44 ± 0.03 


0.900 


9818 


-20.18 ± 0.02 


75.32 ± 2.27 


19.59 ± 0.01 


19.48 ± 0.13 


1.100 


5874 


-20.28 ± 0.03 


69.17 ± 2.90 


19.59 ± 0.02 


19.39 ± 0.14 



Errors on pfit given in the tables are standard errors 
from the fitting procedure, while errors on psuM were calcu- 
lated by dividing the survey areas into sections and finding 
the standard deviation of psum between the sections. For 
the DEEP2 galaxies, we calculated psuM in each of the 4 
fields, whereas for the SDSS uGS galaxies we divided the 
area of the survey into 6 equal areas. The range in mag- 
nitudes psum applies to depends on the magnitude limits 



of the survey and is always less than the fitted magnitude 
range of pfit- As a consequence the fraction psum / Pfit 
ranges from 1.05 to 0.53 for the combined sample of galax- 
ies, 1.05 to 0.63 for the red population, and 1.0 to 0.52 for 
the blue population. This indicates that pfit is never more 
than a factor of 2 extrapolation from psuM- By binning the 
summed LD in 0.2 mag bins and redshift slices, the contri- 
bution of galaxies to the LD, and its evolution as a function 
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Table 4. Schechter function parameters for red galaxies with an assumed a = —0.35. 



< z > 


Ngal 


M* 


4>* (10" 4 gal Mpc" 3 ) 


logp F/T (W Hz" 1 Mpc" 3 ) 


log Psum (W Hz -1 Mpc -3 ) 


0.030 


1708 


-18.60 ± 0.03 


60.91 ±1.74 


18.82 ± 0.01 


18.82 ± 0.12 


0.075 


4414 


-18.69 ± 0.02 


38.80 ± 0.93 


18.66 ± 0.01 


18.65 ± 0.07 


0.125 


4302 


-18.85 ± 0.02 


33.66 ± 1.21 


18.67 ± 0.02 


18.57 ± 0.03 


0.175 


1899 


-19.07 ± 0.03 


25.56 ± 1.71 


18.62 ± 0.03 


18.42 ± 0.08 


0.500 


316 


-19.37 ± 0.08 


20.85 ± 1.62 


18.67 ± 0.03 


18.69 ± 0.08 


0.700 


1296 


-19.44 ± 0.06 


26.06 ± 3.07 


18.78 ± 0.05 


18.71 ± 0.05 


0.900 


1363 


-19.71 ± 0.04 


16.28 ± 0.85 


18.69 ± 0.02 


18.60 ± 0.14 


1.100 


616 


-20.07 ± 0.05 


9.52 ± 0.77 


18.60 ± 0.04 


18.48 ± 0.14 



Table 5. Schechter function parameters for blue galaxies with an assumed a = —1.3. 



< z > 


N G AL 


M* 


(p* (10~ 4 gal Mpc" 3 ) 


\ogp FIT (W Hz -1 Mpc -3 ) 


logpst/M (Wife" 1 Mpc" 3 ) 


0.030 


4466 


-19.03 ± 0.05 


44.50 ± 1.73 


19.01 ± 0.02 


19.01 ± 0.09 


0.075 


9588 


-19.09 ± 0.02 


57.60 ± 1.47 


19.03 ± 0.01 


19.01 ± 0.04 


0.125 


9297 


-19.28 ± 0.02 


53.73 ± 1.86 


19.09 ± 0.02 


18.94 ± 0.03 


0.175 


5941 


-19.51 ± 0.02 


47.85 ± 2.60 


19.12 ± 0.02 


18.87 ± 0.05 


0.500 


1639 


-20.17 ± 0.06 


32.64 ± 1.67 


19.34 ± 0.02 


19.33 ± 0.04 


0.700 


5618 


-20.37 ± 0.04 


33.34 ± 1.24 


19.42 ± 0.02 


19.35 ± 0.03 


0.900 


8455 


-20.39 ± 0.03 


48.23 ± 1.72 


19.59 ± 0.02 


19.42 ± 0.13 


1.100 


5258 


-20.46 ± 0.03 


49.13 ± 2.43 


19.62 ± 0.02 


19.34 ± 0.14 



--X ALL GALAXIES /3 = 1 .36 zt 0.04 
-□ BLUE POPULATION (S = 2.09 4 0.05 
-O RED POPULATION fl3 = -0.26 ± 0,07 
SUMMED RED AND BLUE LD 
BLUE LD ASSUMING « = -' .0 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
Log(l+z) 



Figure 11. Evolution of the M-band luminosity density for the 
different galaxy populations. The LD for the combined population 
is less than the LD for the blue population in the two highest 
redshift slices. This is due to our fixing the slope of the LF for the 
different populations. Black diamonds represent the LD obtained 
from the addition of red and blue LD. Vertical bars indicate la 
Poisson errors on the LD. Horizontal bars represent the redshift 
range. The dashed lines represent fits (eq. to the data. Blue 
crosses indicate the effects of assuming a = —1.0 for the blue 
DEEP2 population. 



of absolute magnitude can be evaluated. Figs. 1 121 1131 and [141 
show this evolution for all, blue and red populations. These 
plots show peaks that indicate we are sampling galaxies that 
contribute the most to the luminosity density. 
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Figure 12. The contribution to the luminosity density for the 
combined population of galaxies for the 8 redshift slices. The LD 
is given in units of 10 18 W Hz -1 Mpc -3 mag -1 . 0.2 mag binning 
is used. 



The plots for the combined and blue populations clearly 
show that as well as the total LD increasing with redshift, 
the galaxies that contribute most to the LD increase in 
brightness with redshift. This illustrates 'downsizing', first 
observed bv lCowie et~all (119961) and also found in more re- 
cent studies ( Treu et al. l l2005l ; lBundv et al.ll2006h . where the 
brighter, more massive galaxies form their stars at higher 
redshift than the lower luminosity galaxies. The LD for the 
red population remains more or less constant over all red- 
shift slices and again the galaxies that contribute most to 
the LD get brighter with increasing redshift. 
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Figure 13. The contribution to the luminosity density for the 
red population. 
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Figure 14. The contribution to the luminosity density for the 
blue population. 



For the combined and blue populatio ns we parametrize 
the evolution in the integrated LD as in iLillv et al l (I1996T I 
by 



p L OC (1 + z) 



(7) 



For the combined population we obtain /3 = 1.36 ±0.04 and 
for the blue population we obtain f3 = 2.09 ±0.05 (assuming 
a constant faint-end slope a). Note if we apply the psum er- 
rors before fitting, we obtain shallower slopes (3 — 1.18 and 
1.95, respectively, and if we apply a = —1.0 to the DEEP2 
blue population (Fig. [TTJ, the evolution is apparently lower 
still. Thus, we estimate there is a 0.2 systematic uncertainty 
in the values of /3 for the combined and blue populations 
(larger for the red population). In the next section we com- 
pare these values of /3 with the results from other surveys at 
different wavelengths. 



5 DISCUSSION 

5.1 Comparison between wavelengths 

First we compare our results with those fro m the /-band se- 
lected VIMOS-VLT Deep Survey (VVDS) (|Le Fevre et all 



|2004 ). Using a sample of 11 034 ga l axies from the first data 
release of the VVDS, lllbert et al.1 <|2005l ) find that M* for 
rest-frame [/-band (360-nm) luminosity function brightens 
by ~ 1.8 mag in the redshift range 0.05 < z < 1. This bright- 
ening is greater than found in our analysis, which shows that 
M* for the c ombined red and bl ue populations brightens by 
~ 1.4 mag. iTresse et al] (|2007h . again using VVDS data, 
find that the evolution in the {/-band luminosity density 
is fitted by (3 ~ 1.9 . By combining the LD measurements 
of lllbert etail (120051) with it-band SDSS main sample data 
from lBlanton et al.l (|2003bh . we find that (3 = 1.6±0.4, which 
is greater than but consistent with our value of j3 = 1.36±0.2 
for the combined galaxy sample. 

In the 150-n m far-ultraviolet (FU V) luminosity func- 
tions produced bv lArnouts et al. (2005) using Galaxy Evolu- 
tion Explorer ( GALEX) satellite data combined with VVDS, 
M* was determi ned to brighten by 2 mag , using a sample 
of 1 039 galaxies. ISchiminovich et al.l l|2005l ). using the same 
data set, find the non-dust corrected l uminosity density to 
evolve with (3 = 2.5 ± 0.7 out to z ~ 1. IWilson et all (|2002h 
found that the 250-nm, UV LD evolves as f3 = 1.7 ± 1.0 
out to z ~ 1.5, using If ,B, V observations taken as part 
of a multi-passband survey of the Hubble Deep Field and 4 
other fields. Our result along with the other more recent UV 
measuremen ts rule out the ear ly estimation of the NUV LD 
evolution bv lLillv et al.l (Il996l ). who found a steep evolution 
of (3 ~ 4 using the Canada- France-Redshift-Survey (CFRS). 

By going to longer wavelengths through optical to the 
NIR, the light sampled in galaxies is emitted by stars of 
increasing lifetimes. For the _B-ban d (440- nm), we deter- 
mined j3 using LD data published in lWillmer et aL (2006). 
This paper made use of an earlier DEEP2 sample of ~ 
11,000 galaxies. By fitting to the data we find the LD 
evolve s as (3 = 0.83 ± 0.46. T his is consistent w i th th e com- 
bined iBlanton et"aH <|2003bh and lllbert et all (|2005h mea- 
surements, which we find gives an evolution in the B-band 
of P = 1.22 ± 0.28. These combined results also allow (3 to 
be estimated for the V, R and /-bands, obtaining j3 values 
of 1.00 ± 0.32, 0.82 ± 0.34 and 0.59 ± 0.34 re spectively. 

At near-infrared (NIR) wavelengths, IPozzetti et al.1 

<|2003t ) found from a sample of 489 galaxies in the range 0.2 < 
z < 1.3, that 13 » 0.7 in the /-band (1.25/im) and (3 w 0.37 in 
the /Cs-band (2.17pim). This /f s -band measurement is also 
consistent with more recent data from the UKIDSS U ltra 
Deep Survey Early Data Release (|Cirasuolo et al.ll2007l ). In 
this paper LFs are produced for 6 redshift slices in the range 
0.25 < z < 2.25 using a sample of ~22 000 galaxies. From 
the LD measured from the first 3 redshift slices, we find that 
(3 = 0.33 ± 0.38 in the range 0.25 < z < 1.25. 

By going to still longer wavelengths in the mid-infrared, 
obscured star formation can be traced from light re-emitted 
by dust grains. From mid-infrared 12 /im data, LFs out to 
z ~ 3 have been produced fro m observations using Spitzer by 
IPerez-Gonzalez et al.1 (|2005l ) for a sample of ~ 8000 galax- 
ies. From these measurements the evolution in LD obtained 
from their own form of the lumin osity function out to z ~ 1 
is best fitted by (3 = 3.83 ± 0.26. iLe Floc'h et al.1 (|2005T ) us- 
ing Spitzer MIPS 24 /im data for 2600 sources, converted 
into total IR luminosities find that the total IR luminosity 
density evolves as (3 = 3.9 ± 0.4. 

Fig. [15] summarises the above results of how the lumi- 
nosity density evolution (/?) varies as a function of wave- 



© 2009 RAS, MNRAS OOO.mfTil 



12 M. Prescott, I. K. Baldry and P. A. James 



length. As wavelength increases from the far UV to the near 
infrared, (3 decreases, which is because of the increasing con- 
tribution from older stellar populations to the luminosity of 
the galaxies. This plot clearly shows the difference between 
our results for the evolution of the it-band LD for the blue 
and combined populations and the evolution of the in frared 
LD (|Perez-Gonzalez et al.ll2005l ; lLe Floc'h et alj|2005l). Also 
shown is the SFR density evolution calculated by Hopkins! 
(2004) using a compilation of X-ray, UV, [O//], H a , H^, 
mid-IR, sub-millimetre and radio measurements corrected 
for dust attenuation where necessary. 



5.2 Estimating the cosmic SFH from the u-band 
LD 

Given the colour bimodality of galaxies, there is a natu- 
ral separation between blue and red populations. We can 
estimate the evolution in the SFR density of the universe 
from the it-band LD evolution by considering just the blue 
population. By doing this we effectively remove most of the 
contribution of it-band luminosity produced by old stellar 
populations (passive red galaxies). 

In order to obtain an estimate of the evolution in the 
SFR density, we first have to consider the residual it-band 
luminosity produced by the old population of stars (r > 1 
Gyr) in blue ga laxies. To do this we perform a n analysis 
using PEGASE (|Fioc fc Rocca-Volmerangelll997l ) models to 
determine the contribution of u-band luminosity from young 
stars wi th ages r < 1 G yr old. In these models we assume an 
IMF of iKroupal (|200lh . with solar metallicity and constant 
star formation rate between a formation redshift (z/ orm ) and 
an observed redshift (z bs)- For a model with zf orm — 6, the 
fraction of the u-band luminosity from young stars changes 
from 0.88 at z b s = 1 to 0.83 at z t s = 0. For zj orrn = 3, the 
fraction changes from 0.9 at z bs = 1 to 0.84 at z b s = 0. The 
assumption that blue galaxies in general form their stars at 
a quasi-con stant rate can be justified considering that stud- 
ies, such as lBrinchmann et alj (|2004T ). have shown low mass 
gala xies to have had a ro ughly constant star formation rate 
and I James et al.l (|2008bh find the same result for local late- 
type galaxies. If the fraction of u-band luminosity produced 
by young stars changes by these amounts from z = 1 to 
z = the evolution in the LD evolution with the contribu- 
tion from old stellar populations removed is (3 — 2.2 ± 0.2. 

The second correction to estimate the evolution in SFR 
density is to consider an increase in dust attenuation from 
z = to z = 1, as luminosities of chara cteristic galaxies at 
higher redshifts are brighter as found bv lCowie et al.l (|l996l ) 
and more luminous/higher SFR galaxies have high er dust 
attenuation than faint systems (|Hopkins et alj|200 ll ). For a 
simple analysis of how much the dust attenuation increases 
from z = to z = 1, we use the magnitudes of the galaxies 
contributing the most to the LD at z — and 2 = 1. Using 
Fig. 1141 it can be seen that the main contribution to the 
luminosity density at z = 1 peaks at M u ~ —20.5 and at 
2 = it peaks at M u ~ —19. Converting these typical galaxy 
magnitudes into flux, we determine the star formation rate 
(in Mq yr" 1 ) of the galaxies , unco rrected for dust, using 
the equatio n of Hopkins ct alj (|200ll ). which is a conversion 
taken from Cra m et al. | il 19981 ) multiplied by a factor of 5.5 
to account for stars with masses ^ 5 Mq- 



SFRuv = L uv I 7.14 x 10 20 WHz" 1 



(8) 



Du st corrected star form ation rates are obtained from 
eq. 7 in iHopkins et al.l l|200l] ) using the SFRs of the typical 
galaxies. The difference between the corrected and uncor- 
rected SFRs gives an effective dust attenuation of A u ~ 1.0 
mag at z = and A u ~ 1.25 at z = 1 which equates to pho- 
ton escape fraction of ~ 0.40 at z = and ~ 0.32 at z = 1. 
This is comparable to the results of lDriver et al.l (|2008l ) who 
find the photon escape fraction to be ~ 0.45 using the local 
galaxies. 

This i s reas onable as other UV studies such as 
iBuat et all (|2005l ) find the mean extinction of a GALEX 
sample of galaxies in the local universe to be A — 1.1 in the 
NUV (230-nm) and A = 1.6 for the 150-nm FUV. The dust 
corrected evolution in the LD for the blue galaxies is essen- 
tially the evolution in the SFR density of the universe ex- 
cluding dusty star forming red galaxies and possibly the most 
extreme objects with hidden star formation such as luminous 
infrared galaxies (LIRGS). Assuming there is an increase in 
the photon escape fraction of 0.32 ± 0.02 to 0.40 ± 0.02 from 
z = 1 to z = 0, this implies that /3s fr = 2.5 ± 0.3. 



5.3 Further correction for dusty star forming 
galaxies 

After correcting for u-band LD evolution for dust and the 
residual u-band luminosity from the old stellar popula- 
tions, the discrepancy still remains between our estimate 
in the SFR evo lution (/3s fr = 2.5 ± 0.3) and the results o f 
IHopkins! (|2004l ) ((3 S fr = 3.3 ± 0.3), iLe Floc'h et~aH (|2005l ) 
{Psfr = 3.9 ± 0.4) and the evolution in SFR densities we 
determined usi n g the combination of UV and far IR from 
Tak euchi et al.l l|2005l ) to give Psfr = 3.4 ± 0.6. Possible 
reasons for the discrepancy include: (i) star formation on 
the red sequence; (ii) a higher dust attenuation correction 
for dusty star forming galaxies; (iii) evolution of the IMF. 

The first way to resolve the discrepancy would be to 
take into account the star formation occ urring in dusty 
galaxi es on the red sequence. Studies such as lFranzetti et all 
(|2007h find, using first epoch VVDS data, that the fraction 
of red galaxies (U — V > 1) that can be spectroscopically 
classified as star forming galaxies with large [On] equivalent 
widths increases with redshift from ~ 10 per cent at z — 0.5 
to ~ 40 at z — 1.1, with 35 per cent of red galaxies being 
spectroscopically classified as star forming overall. Correct- 
ing for these red star-forming galaxies however would only 
change evolution in the SFR by a small amount because the 
contribution to the u-band LD from the red population at 
z = 1 is only of ~ 10 per cent, as seen in Fig. 1111 

Another possibility would be to take into account the 
dust attenuation of starbursting, luminous infrared galaxies 
(LIRGS), fo r which our method for dust correction may be 
inadequate. iPerez-Gonzalez et al.l l|2005l ) find that the con- 
tribution from LIRGS to the SFR density rises from ~ 20 
per cent at z — to become the dominant source of star 
formation at z = 1. In order to obtain an evolution of 
Psfr > 3.2, the photon escape fraction in the u-band would 
have to drop from 0.4 at redshift z = to < 0.2 at z = 1 . 
However, the fossil record analysis of iPanter et alj (|2007l ) 
implies (3sfr ~ 1.5 ± 0.6 (from the first four data points of 
their table 3), which is inconsistent with (3 > 3 (2.5a). 
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Figure 15. Comparison of non dust corrected LD evo- 
lution measured at diff erent wavele n gths. Values from 

ISchiminovich et al.l d2005h. Iwilson et al.l l l2002f) , iBaldrv et all 
l|2005l l and iLe Floc'h et al] 720051) were plotted directly from 
the literatureT TBlanton et all ll2003bl) and lllbert et al.l 11200511 LP 
result s were combined and (5 calculated from a fit. IPozzetti et aD 
l|2003l l is an estimat e quoted in their pa p er and has no errors. 
(} was calcu lated for | Willm er et al] {2006), Perez- Gonzale z" et al.l 
[|2005T ) and ICirasuolo et al .1 j2007l ) using LD values at different 
redshifts and fitting. The error bars for these points are lcr 
errors from the fit. The shaded region indicates the evolution 
of the star form atio n rate estimated from multiwavelength data 
llHopkinsll200i) . For lLe Floc'h et al.l [120051 1 we indicate the range 
of wavelengths that make up the total 1R. 



A final possible way to explain the difference would 
be to consider that starbursting galaxies for m stars with 
a top- heavy IMF, which has been suggested bv lBaugh et alj 
(2005) to explain the number counts of faint sub-millimetre 
galaxies, to explain the metallicit ies of intergalactic gas in 
clusters llNagashima et aD, [2005a) and stars in ellipticals 
l|Nagashima et al.ll2005bh . More recently IL~acev et all (120081 ) 
used a top-heavy IMF to reproduce the evolution in the 
mid-IR luminosity function observed by Sp i tzer. A n evolving 
IMF has been suggested bv lWilkins et al] l|2008l ) to explain 
differences between instantan eous SFRs a nd SFRs inferred 
from the stellar mass density. Dave (2008) also proposes an 
IMF which has more high mass stars compared to low mass 
at earlier epochs in order to explain differences between the 
observed stellar mass-SFR relation and theoretical models. 
It may be possible that the SFR evolution measured from 
the infrared and far UV is more biased to the formation of 
the most massive stars whereas the u-band traces star for- 
mation of a larger mass range of stars in galaxies with a 
'normal' IMF. 



6 CONCLUSIONS 

We have produced and analysed u-band luminosity func- 
tions, to quantify the evolution of the u-band luminosity 
density in the range < z < 1.2, for samples of galax- 
ies taken from the SDSS uGS and DEEP2 galaxy surveys. 



Separating the red and blue populations of galaxies using 
CMDs, our main results are as follows: 

(i) We fit Schechter functions for the red, blue and com- 
bined populations of the galaxies assuming constant faint- 
end slopes of a = —0.35, a = —1.3 and a — —1.0, respec- 
tively. A constant a with redshift gives reasonable fits to 
the blue and combined populations (this does not exclude 
some evolution in a) but is a poor fit to the red sequence at 
z > 0.4, wh ich reflects significant ev olution in the dwarf-to- 
giant ratio (|Gilbank fc BaloghllioOcj ). 

(ii) M* has brightened by ~ 1.4 mags for the combined 
popula tion from z = to 1.2, which is similar to the find- 
ings of lllbert et al.l (|2005l ) using the VVDS sample who find 
galaxies brighten by 1.6 to 2 mags. 

(iii) The galaxies in all populations that contribute the 
most to the M-band luminosity are observed to increase in 
bri ghtness with redshi ft. This illustrates 'downsizing' as seen 
bv lCowie et alj (|l996l ) and later studies. 

(iv) By parametrizing the evolution in u-band luminosity 
density as p cx (l+z) 13 , we find that the combined population 
of galaxies evolves with f3 — 1.36 ± 0.2 and for the blue 
population f3 = 2.09 ± 0.2. The red population LD remains 
approximately constant over time. 

(v) Considering just the blue population, removal of the 
u-band luminosity contribution from old stellar populations, 
estimated from population synthesis models, increases the 
u-band LD evolution to /3 — 2.2. This represents the SFR 
density evolution uncorrected for dust. 

(vi) We estimate that the average dust attenuation is 1.0 
mag at z = and 1.25 mag at z = 1. By correcting the 
u-band LD for this extinction, we obtain the evolution in 
the SFR excluding red star-forming galaxies and LIRGS to 
be /3sfr = 2.5 ± 0.3. This modest correction for evolution 
in dust attenuation may be appropriate for estimating the 
build-up of stellar mass, whereas the more severe dust evo- 
lution suggested by mid and far-IR measurements may be 
biased toward top-heavy IMF star formation. 
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